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ABSTRACT 

2 

For  a  sequence  Z  :=  { z . ..,  Z|z|)  in  Z  ,  the  h°x  spline  Mz  is  the 
distribution  given  by  the  rule 

MZ$  :=  /  K5*  d*  * 


:  i' 


We  show  that  P„(x)  :=  Z  „  M(k)  elkX  >  0  iff  AZ  C  {(1,0),  (0,1),  (1,1)}  for 

2  Z  — 

kes 

some  linear  map  A  which  is  1-1  on  Z2.  This  implies  that,  for  such  a  Z  ,  there 

exists,  for  any  bounded  continuous  function  f  ,  a  unique  bounded  spline 

Izf  e  span{Mz(*-k)  s  ke*2}  which  agrees  with  f  on  Z2  .  We  give  sufficient 

conditions  in  terms  of  the  Fourier  transform  of  f  for  the  convergence  of  Izf  as  the 

degree  | Z | — 2  tends  to  infinity.  E.g.,  we  show,  for  the  case  tnat  Z  contains  each 

of  the  vectors  (1,0),  (0,1),  (1,1)  exactly  n  times,  that  If  -  I  fl  — TTL>  0 

n  -,R  " 

A  9  ** 

supp 


f  CC  fl  conv{ (z,z) ,  -(z,z),  (-z,2z),  -(-z,2z),  (-2z,z),  -(-2z,z)}  ,  with 


2n 

z  :=  —  .  The  converse  of  this  result  holds  with  "CC"  replaced  by  "C" 

J  y 
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SIGNIFICANCE  AND  EXPLANATION 


'' in  a  series  of  monographs^ (cf.  MRC  Technical  Sunmary  Reports  #852,  #1104, 
#1113,  #1150,  #1153,  #1180, '#1184,  #1219,  #1222,  #1326,  #1330,  #1360,  #1423, 
#1546)  I.  J.  Schoenberg  developed  a  comprehensive  theory  of  univariate 
cardinal  splines.  His  results  strongly  influenced  the  analysis  of  totally 


positive  matrices. 

7, 

In  this  report  wT extend  two  of  his  basic  results  on  cardinal 

interpolation  to  bivariate  box-splines, which  have  been  introduced  in  MRC 

'  -r*  j  i 

Technical  Summary  Reports  #2320,  #2415.  show  that,  for  functions  of 
exponential  type,  cardinal  interpolation  is  a  rapidly  convergent  approximation 
process  as  the  degree  tends  to  infinity.  Being  not  restricted  to  a  tensor 
product  mesh  gives  a  greater  flexibility,  and  because  of  the  exponential  decay 
of  the  Lagrange  functions,  spline  interpolation  is  suitable,  e.g.,  for  data 

.—r 

A-  y 

smoothing.  WB'  also  expect  that  bivariate  cardinal  splines  have  a  similar 
significance  for  theoretical  questions  as  in  the  univariate  case.  \ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumnary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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to  whose  insight  and  sense  of  beauty  we  are  all  indebted 


1.  Introduction,  in  this  paper,  we  carry  Schoenberg's  beautiful  cardinal  spline 
theory  [Sj  3]  over  to  a  two-dimensional  context  which  is  not  just  the  tensor  product  of  the 
univariate  situation.  We  find  that  we  must  work  harder,  yet  must  be  satisfied  with  less 
precise  results. 

He  are  after  a  bounded  cardinal  interpolant  to  bounded  data.  This  means  that  we  are 
looking  for  a  function  of  the  form 


If 


Z  a  M( «-j  ) 

je*  3 


with  a  e  which  agrees  with  a  given  bounded  function  f  on  I2  .  Here,  M  is  a 

fixed  function  of  compact  support.  In  Section  2,  we  follow  Schoenberg  [S3 ]  in  describing 
necessary  and  sufficient  conditions  on  the  Fourier  transform  of  M  to  insure  the 
correctness  of  the  interpolation  problem,  i.e.,  the  existence  and  uniqueness  of  solutions. 

We  are  particularly  interested  in  using  for  M  a  box  spline  ,  i.e.,  the  two- 
dimensional  "shadow"  of  an  m-dimensional  cube,  as  given  explicitly  in  (1.1)  below.  We  find 
it  convenient  to  change  the  definition  [ BH 3 ] 


M  *  :=*  /  ,  ♦(  S  AUK) 

(o,n  cez 

of  the  box  spline  M  »  Mz  to  include  an  appropriate  shift  which  makes  the  origin  the 
center  of  the  support  of  M  .  This  means  that  we  use  the  definition 
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M« 

This  gives  the  Fourier  transform 


with 


s=  /  ♦(  E  X(C)C)  dA  . 

(J/2  ,V2  )  cez 

M*  of  M  the  symmetric  form 

M*(x)  -  n  S(C*x) 

cez 


(1.1) 


(1.2) 


Z 


S(t) 

It  is  obvious  from  this  formula  that 
are  replaced  by  their  negative,  i.e.. 


sin  t/2 
t/2  * 

M  =  Mz  is  unchanged  if  one  or  more  of  the 


0.3) 

c  e 


maz  -  H*  H-4) 

if  A  -  dlag(±1, . . . ,±1 )  .  Further,  if  A  is  any  matrix,  then 

M*  (x)  =■  M*(A*x)  and  M  (Ax)  -  M  (x)/det  A  .  (1.5) 

AZ  Z  AZ  Z 

This  allows  one  to  deduce  symmetries  in  M*  in  case  AZ  equals  Z  after,  possibly,  some 
elements  of  AZ  have  been  multiplied  by  -1  . 

The  set  Z  of  directions  can,  of  course,  be  chosen  arbitrarily.  But  since  we  are 
interested  in  having 

S  :»  span  (w(*-j))  ? 

je* 

be  a  simple  piecewise  polynomial  space,  we  choose  Z  from  X2  .  It  is  shown  in  [BH1 ]  that 
the  integer  translates  HI*-])  ,  j  «  I2  of  the  box  spline  are  linearly  dependent  (when 
allowing  for  infinite  linear  combinations)  in  case  the  direction  set  Z  contains  two 
vectors  which  span  a  proper  sublattice  of  Z2  •  Linear  independence  is  an  obvious  necessary 
condition  for  the  cardinal  interpolation  problem  to  be  correct.  Thus,  up  to  obvious 
symmetries,  this  leaves  the  three  vectors  (1,0),  (0,1)  and  (1,1)  as  the  only  candidates 
for  the  directions  C  in  Z  . 

With  this  restriction,  S  is  a  space  of  piecewise  polynomial  functions,  of  polynomial 
degree  I Z I —2  or  less,  and  with  possible  discontinuities  only  across  the  three  types 
x(1)  -  k  ,  x( 2 )  -  k  ,  x(1)  -  x(2>  -  k  ,  k  e  * 


of  mesh  lines.  The  overall  smoothness  of  the  elements  of  S  depends  on  the  multiplicities 
of  the  directions  in  Z  .  Such  details,  as  well  as  the  relationship  of  S  to  the  space  of 
all  piecewise  polynomial  functions  on  such  a  three -direct ion  mesh,  of  degree  I Z | —2  and  of 
specified  smoothness,  are  all  discussed  in  [BHp) • 


-  1.2 


In  Section  3,  we  supply  certain  detail  concerning  symmetries  of  such  a  three-direction 
box  spline  and  its  Fourier  transform,  we  prove  the  correctness  of  cardinal  interpolation 
with  such  a  box  spline  in  Section  4.  We  spend  the  major  effort  of  this  paper  in  Section  5 
where  we  prove  that,  under  reasonable  conditions,  the  cardinal  interpolant  If  of  any 
suitably  smooth  function  f  converges  to  f  as  |z|  — >  •  .  Specifically,  we  prove  such 
convergence  under  the  condition  that  f  is  the  Fourier  transform  of  some  compactly 
supported  measure,  following  entirely  the  path  established  by  Schoenberg  [S]  in  the 
univariate  case  who  showed  that  such  convergence  could  be  had  whenever  supp  f *  c  (-»,»)  . 
We  find,  though,  that,  in  our  bivariate  setup,  there  are  many  different  sets  playing  the 
role  of  this  interval,  and  which  of  these  sets  is  relevant  depends  on  the  manner  in  which 
I Z |  goes  to  infinity. 

The  final  section  is  devoted  to  the  many  detailed  estimates  on  which  the  arguments  in 


Section  S  are  based, 


2.  Cardinal  Interpolation.  Let  M:R2  — >  R  be  a  continuous  function  with  compact 
support,  and  denote  by 

S  :«  Sjj  :=»  span  {M(  •  -  j)  s  j  6  Z2} 

the  space  generated  by  its  integer  translates.  Cardinal  interpolation  with  M  concerns 
inversion  of  the  linear  map 

5  'y  ■ 

We  say  that  cardinal  interpolation  with  M  is  correct  if  this  map  is  1-1  and  onto,  hence 
boundedly  invertible,  and  denote  its  inverse  by  IM  or  I  .  In  other  words,  cardinal 

2 

Interpolation  with  M  is  correct  iff  there  exists,  for  every  bounded  sequence  f  e  1^(2  ) 
,  a  bounded  function  If  e  S  which  agrees  with  f  on  Z2  .  The  interpolation  problem, 
i.e.,  the  determination  of  If  ,  is  equivalent  to  the  algebraic  problem  of  determining  the 
coefficient  sequence  a  for 


so  that 


If  -  E  a^MC- J) 


2 

a  e  1^  and  E  a^M(»-j)  ~  f  on  Z  . 

Hence  the  correctness  of  cardinal  interpolation  is  equivalent  to  the  lnvertibility  of  the 
matrix 


A  s-  («<  3->'>)jfkex2  (2.2) 
as  a  map  on  1—  .  Since  A  is  a  banded  (bivariate)  Toeplitz  matrix,  we  have  the  following 
necessary  and  sufficient  condition  for  the  correctness  of  cardinal  Interpolation. 


theorem  2.  Cardinal  interpolation  with  M  is  correct  iff 


P(x) 


PM(x)  j-  E  M( j )  e 


ijx 


(2.3) 


does  not  vanish. 


Proof.  If  P(x)  “  0  ,  then  (e  e  ker  AH  i  ,  and  this  contradicts  the 

assumption  that  A  is  1-1.  On  the  other  hand,  if  P  does  not  vanish,  then  the  inverse 
of  A  can  be  expressed  as  a  Toeplitz  matrix. 


(A-1) 


jk  " 


/ 


<-!»,*) 


-i( j-k)x 
P(x) 


dx/2e  . 


(2.4) 


2.1 


(2.5) 


In  view  of  the  geometric  decay  of  the  Fourier  series  for  1/P  ,  we  have 

|(A-1)jkl  <  const 

for  some  A  »  A(P)  e  (0,1)  .  Therefore,  A-1  is  bounded  on  l  (Z2)  for  any  p  @ 

P 

[1,-1  -III 


It  is  convenient  to  write  the  cardinal  interpolant  in  Lagrange  form: 

If  =  If.  L<—  j) 

1 


L  16  -  Z  (A~  )Q..  M(»-j)  (2.6) 
the  finnl— ni si  function  of  the  interpolation  process.  The  Fourier  transform  L*  of  L  is 
particularly  simple.  Combining  (4)  with  (6),  we  obtain 

L*  -  M*/P  .  (2.7) 


He  will  also  make  use  of  the  identity 


P(x)  *  Z  M*(2»j  -  x) 

which  follows  from  applying  the  Poisson  summation  formula  Ef(j)  =  Ef“(2wj)  to  (3). 


3.  Cardinal  interpolation  with  a  box  spline.  In  this  section,  we  develop  in  some 
detail  facts  about  cardinal  interpolation  with  the  box  spline  M2  .  Recall  from  Section  1 
that  (m(#  —  j ))  is  linearly  dependent  if  Z  contains  two  vectors  which  span  a  proper 

je*r 

sublattice  of  TT  .  Linear  independence  of  (m( •— j ) )  is  an  obvious  necessary  condition  for 
cardinal  interpolation  with  M2  to  be  correct.  Thus,  up  to  obvious  symmetries,  the  only 
relevant  case  to  consider  is  the  case  when  the  only  directions  in  Z  are 
dj  :=  (1,0),  d2  :=  (0,1)  ,  and  d3  :=  (1,1)  . 

We  show  in  Section  4  that,  with  this  restriction,  cardinal  interpolation  with  Mz  is 
always  correct. 

Assume  from  now  on  that 


Z  -  (d1:r,  d2:s,  d-jit)  . 

In  this  case,  Z  is  characterized  by  the  vector 

n  :=  (n1#  n2,  n3)  (r,s,t) 

of  direction  multiplicities,  and  we  will  freely  write  n  instead  of  Z  whenever  it  is 

necessary  to  indicate  by  subscript  the  dependence  on  Z  of  some  quantity.  Further,  the 

general  formulae  given  in  Section  1  simplify.  For  example, 

M*(u,v)  »  S(u)r  S(v)s  Slu+v)*  ,  (3.1) 

n 

with 

s<t>  .  t/2 

Further,  the  characteristic  polynomial  P  *  Pn  and  the  Fourier  transform  L“  of  the 


fundamental  spline  L  *  Ln  have  the  representations 

>1 


P(2wu,2sv)  =■  *,n|  (sin(xu))r(sin(xv))B(sin(ir(u+v)))t  l  - — 


rk+sl+t(k+t) 


(u+k)‘  (v+l)S(u+v+k+t)t 


(3.2) 


and 


i/L*  (2xu,2sv)  -  E  (-)rk+st+t(k+t,(^:)r(-T7n-— * 

vu+k '  '•v+l'  '■u+v+k+t' 


#  k,tes 

Let  A  m  denote  the  transpose  of  A  .The  relation 


(3.3) 


MAZ<x)  ‘  MZ(AX) 

valid  for  any  matrix  A  together  with  the  fact  that 


(1.5) 


maz  =  Mz 


(1.4) 


in  case  A  -  diag(e  , . . .  ,e  .  .)  with  e,  e  {-1,1} 
i  I  n  |  i 


all  i  ,  implies  certain  syimnetries  of 


M  and  M*  if  the  matrix  A  leaves  the  set 


Z±  !=  {<11*d2'd3'  -d1*_d2'-d3^ 

invariant.  Denote  by  A  the  group  of  all  such  invertible  matrices  A  .  Each  A  €  A  is 

associated  with  a  permutation  oft  e  $3  ( :=  symmetric  group  on  3  elements)  by  the  condition 

Adi  e  {doACi )'  -\(i>}  '  i’1'2'3  * 

From  the  two  matrices  corresponding  to  a  given  06)31  we  choose  one,  A^  ,  in  such  a  way 

that  the  six  matrices  form  a  group  and  we  call  this  group  A+  .  Thus, 

Vi  *  {d0(i)'  -dO(i)}  '  “U  °e$3  '  (3-4) 

and  one  choice  for  the  group  generators  are  the  three  matrices 


f  0 

1-1  o'  ' 


‘(13) 


(  ;  .?)  - 


•=  r1  ’i 

1 0 


corresponding  to  the  transpositions  (12),  (13),  and  (23)  .  With  the  definition 

o(n)  s-  <na(1)»  n0(2)»  "bO)’  ' 
it  follows  from  (1.S)  and  from  (4)  that 

Vx)  -  Mo(n)<±V)  ' 

Mi(n,(v>  -  • 

This  implies 

V^V’  “  P0(n)(x> 

and 


(3.5) 


(3.6) 


1 


I 


u 


■s 


:=  {A  :  A  e  A} 


in  more  detail.  Set 


-•  -  (.!  J) d ' 

Z't  {d’:  d  e  Z±}  -  {(0,1),  (-1,0),  (-1,1),  (0,-1),  (1,0),  (1,-1)}  . 

Since  d**d  «  0  ,  we  see  from  (4)  that  A*  leaves  Z^  invariant.  To  further  illustrate 
the  action  of  the  group  A*  ,  we  divide  f?  into  the  six  cones  Rg  ,  o  e  ,  as 

indicated  in  Figure  3.1. 


Figure  3.1. 

It  is  easily  checked  that 

* 

A  R  *  R  ,  all  o  e  $,  .  (3.8) 

0  0  3 


-  3.3  - 


4.  The  correctness  of  cardinal  interpolation  with  .  In  this  section,  we  show  that 
cardinal  interpolation  with  Mn  is  correct  for  all  choices  of  n  e  S+3  . 


Theorem  4.  For  all  n  e  8+3  ,  Pn  is  strictly  positive. 


Since  P  is  2r -periodic,  this  amounts  to  the  claim  that 

Pn(x)  >  0  for  all  x  e  [-»,*]  2  . 

This  is  the  bivariate  analogue  of  Schoenberg's  well  known  result  for  univariate  cardinal 
spline  interpolation.  To  recall  this  result,  denote  by 


the  univariate  cardinal  B-spline  of  degree  r  ,  and  by  Qr  the  corresponding 
characteristic  polynomial  given  by 

ijx 


Qr(x)  s-  I  Nr(j)  e 


Schoenberg  showed  in  [Sj]  that 


mlnx  2r<x>  “  2r(1,)  “  2(»)r+1 


(-) 


v(r+1 ) 


(4.1 


(2v+1) 

The  fact  that,  for  any  r  ,  the  minimum  is  attained  at  x  «  x  is  a  consequence  of  the 
total  positivity  of  the  matrix  (n  (j-k)).  .  .  In  view  of  this  result,  one  might  think 

r  j  t KcZ 

that,  in  the  above  theorem,  rain^y  P(u,v)  ”  P(w,ir)  .  This  is  trivially  true  in  the  tensor 
product  case,  i.e.,  when  n  «  (r,r,0)  .  However,  in  general,  the  point  at  which  Pn 
attains  its  minimum  depends  on  n  .  It  would  be  interesting  to  determine  its  location  for 
special  choices  of  n  .  The  nicest  conjecture  in  this  context  (cf.  Section  5)  is  that 

min  P(u,v)  ”  P(2x/3,  2w/3 )  in  case  n  «  (r,r,r)  .  (4.2 

In  the  proof  of  the  theorem,  we  make  use  of  (3.6).  This  allows  us  to  assume  without 
loss  of  generality  that 


r  >  s  >  t  . 

We  first  consider  two  cases  which  reduce  to  Schoenberg's  result. 

The  tensor  product  case  n  »  (r,s,0)  .  Here,  we  have  Mn(u,v)  •  Nr(u)Na(v)  ,  and  this 
implies  that 


Pn(u,v) 


Qr(u)Q,(v) 


The  case  n  »  (r , 1 , 1 )  .  Since  the  open  support  of  Mn  intersects  exactly  one  mesh 

line  of  the  form  (*,4)  ,  viz.  the  raeshline  (*,0)  ,  it  follows  that,  in  this  case, 

IN  (k)  .  t  -  0 
r 

o  ,  lei  \  o 

This  means  that  cardinal  interpolation  with  Hn  reduces  to  univariate  interpolation  with 
Nr  on  each  of  the  lines  (*,t)  ,  t  S  I  .  In  particular,  P(u,v)  -  Qr(u)  . 

For  the  proof  of  Theorem  3,  it  remains  to  consider  the  cases  where  the  multiplicities 
are  all  at  least  1  ,  with  equality  for  at  most  one.  We  make  this  assumption  for  the 
remainder  of  this  section. 

To  prove  the  positivity  of  P  ,  we  use  the  representation  (2.10)  in  the  form 

P(2*x)  -  I  K~(x+j)  ,  (4.3) 

with 

M~(x)  M*(2fx)  . 

Recall  from  (3.1)  that,  for  x  »  (u,v)  and  j  «  (k,t)  , 

,  ,  .  ,  .rk+s*+t(k+Jl) 

M  (x+j)  «  w  n  (sin  »u)r(sin  »v)S(sin  »(u+v))  - —  .  (4.4) 

(u+k)r(v+f)E(u+v+k+t)t 

It  is  sufficient  to  show  the  positivity  of  P(2»*)  on  [0 , 1/21 2  for  arbitrary  n  . 

This  follows  from  (3.6)  since,  by  (3.8), 

[-1/2, 1/2] 2  C  LI  A* (0 , 1/2] 2  . 

~  Ae* 

For  x  e  [0,1/2]2  ,  we  now  show  that  the  three  positive  terms 

M  (x)  ,  M  (x  -  d^  ,  and  M  (x  -  d^)  (4.5) 

dominate  the  sum  in  (3).  To  this  end,  we  associate  each  of  the  other  terms  with  one  of 
these  (even  to  the  point  of  splitting  one  of  the  other  terms  between  two  of  these)  and  show 
that  the  resulting  three  sums,  when  divided  by  their  respective  dominant  term,  each  is  less 
than  1  .  For  ease  of  argument,  we  actually  split  the  sum  into  altogether  ten  parts,  as 


indicated  in  part  by  the  figure  below 


Figure  4. 1 


To  simplify  notation,  we  set 

b  (j>  b„  <j.x)  s-  |M  -X^>l  ,  v-1,2,3  ,  (4.6 

W  W'"  M~( x+ j  ) 

with 


Jl  !“  0  »  j2  s*  -<J1  •  33  ~d2  • 

We  now  prove  that 


EjeJ,uJ4 

b,(l>  +  s*#o  bi 

u. 

-t)  < 

1 

(4.7 

jea2uj5 

b2(J)  ♦ 

EU|>o  V'1'1’ 

♦ 

b2(-1. 

-1)/2  <  t 

(4.B 

jej3uj6 

b3(j)  ♦ 

EUI>,  V1'"” 

■f 

b3(-1. 

-1)/2  <  1 

(4.9 

Since  each  of  the  summands  (divided  by  its  appropriate  dominant  term)  other  than  the  three 


dominant  terms  (5)  occurs  in  (7-9)  exactly  once,  we  conclude  from  (7-9)  the  positivity  of 


The  estimation  of  the  various  sums  in  (7-9)  is  straightforward.  In  each  case,  we  find 


a  majorant  which  is  independent  of  x  e  [0,1/2] 2  and  n  .  For  this,  recall  that  we  are 

assuming  that  r,s,t  >  1  with  at  most  one  equality. 

We  begin  with  the  sum  Z  b  (j)  .  By  (3. 1,4. 6),  we  have  for  x  =  (u,v)  and  j 

J1  1 


(k,t)  , 


r  s  t 

.  ,  j.  i  u  ,  ,  v  ,  |  u+v  | 

b1  J  =  u+k  v+ 1  u+v+k+t  ' 


Since  k,  t  >  0  for  j  =  (k,l)  e  and  we  are  assuming  that  u,  v  e  [0,1/2]  ,  this 
quotient  is  largest  when  u  «  v  »  1/2  ,  i.e., 

b1(  j)  <  ©Wfi)1' 

This  bound  is  largest  when  the  exponents  r  ,  s  ,  t  are  as  small  as  possible  ,  i.e., 
when  n  =  (1,2,2),  (2,1,2)  or  (2,2,1)  .  Since 


E  r  v2  i  r  v2-w  r_i  .  1723 

k , t>0  M/2+k-1  1/2+t'  h+k+l^  *  1723.  . 

(k,i)*0 


for  these  values  of  n  ,  we  conclude  that 


b,(j>  <  .18 

J1  1 


Similarly,  one  verifies  that,  for  j  =  (-k,-£)  e  J4  and  x  €  [0,1/2]*  , 

b,d)  < 


(4.10) 


and  so  obtains  that 


E.  b  ( j)  <  .02  , 

J4  ' 


Z  (2k- 1 )”r ( 21- 1 )”8 (k+l-1 ) _t 


.0101. .  . 


k,t>1 

for  (r,s,t)  -  (1,2,2),  (2,1,2),  or  (2,2,1). 


Finally,  for  j  ■  (-t,t)  and  l  *  0  ,  we  have 


b^-t.t)  < 


f-J/2-)1  f-IZi-]8  t  >  o 

^t+1/2 '  U-1/2j 


and  so  obtain 


(4.12) 


Combining  (10-12)  establishes  (7). 

The  other  inequalities  are  proved  in  a  similar  fashion  and  we  only  list  the  estimates 
involved. 


Proof  of  (8): 


b2(j) 


,  r  s  .  ,  t 

Hlli  i_v_|  ,_u+Yll_| 
u+k'  ‘v+t  'u+v+k+t' 


j  -  (-k,t)  e  J2 


b2(j) 


< 


1/2 

t+1/2 


and 


Therefore 


m  lr—  1 

Ek=2£t=0 


1/2 

t+1/2 


.25... 

.23...  for  (r,s,t)  = 
.  30 . . . 


(2,2,1) 
(2,1,2)  . 
(1,2,2) 


j  -  (k,-l)  e  Jc  : 


ET  b  ( j )  <  .35  . 

J2  2 


V”  < 


and 


Therefore 


s  .  t 


.027. .. 


Wt«2  (£)  (Fl72^  (*=*)  “  for 


.010. 


(2,2,1) 

(2,1,2) 

(1,2,2) 


j  «  (-1,t)  ,  |t  |  >  1  s 


b2(j)  < 


IT  b,(j)  <  .1 

J5  2 


(^)8(rZJr)t  .  »<-i 

V2  w  1  >2  .  f  1/2  W  1 


Ithl  b2<j>  <  rt>1^t  +  1/2^t- J  +  U-1/2^t+J  1 


.329.  . . 


j  -  (-1,-1): 


b2(-1.-1>  <  71751^  <  1/3  (at  (u,v)  *  (0,1/2)). 


Proof  of  (9)s  Since 


M?  =  k|(u'V>  =  .  .  ( V,  u ) 

(r,s,t)  (s,r,t) 


and  therefore 


b2,(r,s,t)((u'v,'(k,1))  “  b3, (s,r,t)((v'u)'(t'k))  ' 


4.5 


the  inequality  (9)  follows  from  (8)  by  interchanging  the  roles  of  u  and  v  as  well  as 
those  of  k  and  t  . 

This  completes  the  proof  of  Theorem  3.  |  |  | 


5.  Oonnrgvnc*  of  Cardinal  Interpolation.  'Ms  section  is  devoted  to  the  main  goal  of 


our  paper,  a  study  of  the  convergence  of  the  cardinal  interpolant  to  smooth  functions  as 
the  degree  tends  to  infinity.  We  prove  the  analogue  of  I.J.  Schoenberg's  basic  result: 


Theorem  [82]*  .If  f  is  the  Fourier  transform  of  a  measure  with  support  in  (-*,»)  , 
then  Its  cardinal  spline  Interpolant  Irf  of  degree  r  converges  to  f  as  the  degree 
tends  to  infinity,  i.e.. 

If  -  I  fl  - >  0  . 

r  “  r*“ 


This  theorem  is  a  consequence  of  the  fact  that  the  Fourier  transform  L*  of  the 
fundamental  spline  converges  to  the  characteristic  function  of  the  interval  (-«,*) 


The  bivariate  situation  is  more  complicated.  Here,  the  limit  of  L“  depends  on  just 

n 


how  n  goes  to  infinity.  Recall  from  (3.3)  that 


1/L*(2*x) 


1  ♦  E  e  (x)  a  (x)  , 
je*2\o  3  n'3 


(5.1) 


with  e  (x)  e  {-1,1}  and 

^  r  b  t 

/  %  l  u  I  I  v  .  ,  u+v  , 

*n,  (k,t)  ^U'V  «“  W  W 


(5.2) 


Define  the  "middle  coogjonent"  |i(n)  of  n  by  the  requirement  that  it  equal  the  middle  or 
second  number  in  any  ordering  of  the  three  numbers  r,  s,  t  ,  and  set 

n'  :»  (r'.s'.t')  :=*  n/li(n)  . 


Then,  the  typical  summand  in  the  right  hand  side  of  (1)  is,  up  to  sign, 

>M(n) 


an,j(x)  ‘  (V,j(x)) 

This  shows  that  L*(2«x)  is  close  to  1  for  large  u(n)  provided  an<  j(x)  <  1  for 

all  J  *  0  . 


The  set 

{x  :  a  ,  (x)  <  1  for  all  j  e  Z2  No)  . 

n  ,  j 

depends  on  n'  .  In  particular,  we  cannot  expect  it  to  converge  as  |n|  — >  *  unless  n'  » 
n/u(n)  converges,  to  some  3-vector  m  ,  say.  Here,  we  are  willing  to  allow  m  to  have 

2 

infinite  components.  For  example,  if  n  =  (1,s,s  )  ,  then  p(n)  ■  s  jnd  n*  = 

(1/s, 1,s)  — >  (0,1, “)  as  s  — >  •  .  But  not  every  m  e  [0,«13  is  a  possible  limit.  By 


J  Z'±  -  {(1,0), (0,1). (-1,1), (-1.0), (0,-1), (1,-1)} 

play  a  major  role.  Note  that 

Q  -  0  ,  . 

n  n* 

A  qualitatively  correct  picture  of  is  given  in  Figure  5.1  which  shows  the  roughly 

hexagonal  shape  of  fl  and  also  shows  the  six  curves 
n 

r  {x  e  C  !  a  (x)  =  1}  ,  j  e  J 

n.J  “3  n,j 

which  contribute  to  the  boundary,  3Q  .  Here,  C  4  is  the  union  of  the  two  cones  R 

n  j  o 

which  contain  -j  . 

Themes  5«1»  For  m  e  N  ,  let  be  the  characteristic  function  of  .  Then,  for 

any  d  >  0  ,  there  exists  e  >  0  so  that 

|L'(2*x)  -  x_(x)|  <  C.fl  +  C  diat(x,30  (5.3) 

n  m  1  m 

for  all  x  with  dlst(x,  id  )  >  d  and  all  n  e  (0,«)3  with  n'  t”  n/u(n)  e  B  (m)  ,  and 
■  "r“  m  ■■■"“  c 

with  the  positive  constants  and  C  Independent  of  m  ,  n,  d,  or_  x  . 

Proof.  The  proof  is  baaed  on  a  aeries  of  propositions  which  we  merely  state  as  needed 
and  prove  at  leisure  later.  Ms  begin  with  the  following 

Propositi am  5.1.  fln  depends  continuously  on  n  in  the  Hausdorff  topology. 

which  is  part  of  the  Corollary  to  Lemma  6.4  below.  This  implies  ,  given  d  >  0  ,  the 
existence  of  e  >  0  so  that 

dist(n  ,,n  >  <  d/2 

n'  m 

for  all  n'  e  B£(ra)  .  Consequently, 

dist(x,3H  )  <  2dist(x,3n  ,) 

in  n 

for  all  x  with  dist(x,9fl  )  >  d  and  for  all  n'  e  B  (m)  .  It  is  therefore  sufficient  to 

n  c 

prove  (3)  with  m  replaced  by  n'  . 

For  its  proof,  we  use  (1)  and  we  consider  two  cases. 

(i)  x  e  .  He  need 


-  5.3 


Proposition  5.2.  Lot 


J’  i-  A  (1,1)  -  {1(1,1),  ±(2,-1),  ±(-1,2)} 


Fog  n  8  M  and  x  e  Qn,  “  Qn  , 


(l  +  C  dlat(x,  30  ))"  ,  j  e  JUJ 

'J  (1  +  C  131]"  ,  J  e  a2  \ 


,  j  e  a2  \  (o ujuj1  ) 


with  the  positive  conatant  C  Independent  of  n  ,  j  ,  or  x 


This,  together  with  (1),  implies  that 


11/t;  (2SX)  -  1|  <  Ej#0(«n.<j(x)),‘ 


<  12(l  +  C  diat(x,  30  1))"‘,(n’  +  I  (1  +  C|j|)"M<n)  (5. 

n  jexoujur' 


<  ct(i  +  c  disttx,  3fln, )]' 
and  so  proves  (3)  for  this  case. 

(11)  x  f  9  ,  •  For  this  case,  we  need 
n 


Proposition  5.3.  The  Integer  translates  of  0  fora,  up  to  a  set  of  measure  zero,  a 


tltlon  of  A2  ,  l.e. , 


K2  -  U,  3  +  al  •  anu  +  a)-fi  for  J*o. 
3«2  n  n  n  - 


Me  conclude  that  there  is  )  /  0  so  that  x  ■  x’  +  3  with  x'  e  .  With  this, 
the  periodicity  of  the  characteristic  polynomial  p  to  write  for  such  x 


L*(2wx)  -  L* (2» (x'+3 ) ) 


M  (x*+3)  M  (x')  _  ..  ..  M  (x'+3  > 

P(2w(x'+1))  -  1  '  ~ 

'  M  (X1)  M  (X*) 


Therefore 


By  (5), 


|L*(2»x)|  -  |L*(2»x’)l  a  .(x’)  -  |L*(2wx’ ) | (a  ,  .(x'))M,n)  . 

n  n  n,3  n  x  n  ,3 


|L*(2wx* ) |  <  C 

n 


Thus  (3)  is  proved  for  this  case  once  we  show 


Proposition  5.4.  Let  x  «  x*  +  j  with  j  e  ^  \  0  and  x'  e  ,  and  with  n  e  R3  . 

Then 

an,(j(x')  <  {l  +  C  dist(x,aQn, >)-1  (5.6) 

for  some  positive  conatant  C  Independent  of  n  and  x  . 

This  finishes  the  proof  of  Theorem  5.1.  ||| 

Theorem  5.2.  Let  f  he  the  Fourier  transform  of  a  measure  with  support  strictly 
inside  2xfl  for  some  m  e  N  ,  i.e., 

d  diat( supp  f",3(2*fl  )  )  >  0  . 

IB 

Then  there  exists  e  >  0  so  that,  for  all  n  e  B£(m)  , 

If  -  I  fl^  <  Ct(l  ♦  Cd)"M<n) If *I1  . 

with  ****■)  the  total  variation  of  f*  .  ’CTie  positive  constants  C  ,  C1  do  not  depend 
on  m  ,  d  ,  or  n  . 


Proof.  Fix  d  >  0  and  choose  C  >  0  so  that  supp  f *  C  (2*0^,)  And 

diet (supp  f * ,  3(2*0^,))  >  d/2  for  all  n*  6  B£(m)  .  He  have  to  estimate 

f(x)  -  I  f(j)  L  (x-j)  -  f<x)  -  r(2»)"2  /,  f( j)e'ljyL*(y)  elxydy  . 

n  n 

Since  2iBn,  is  a  fundamental  domain,  i.e..  its  integer  translates  form  a  partition  of 

unity  (by  Proposition  5.3),  and  supp  f*  C  Q  ,  ,  (f(j))  are  the  Fourier  coefficients 

of  the  periodic  extension  f  of  the  measure  f “  .  Using  the  weak  convergence  of  the 

Fourier  series  of  a  measure,  we  obtain 

f(x)  -  (X  f)(x)  -  <2*)"2/,  [f “(y)  -  f*  (y)  L*(y)l  eixy  dy  . 

n  *•  p  n 

Applying  Theorem  5.1  yields,  for  n'  e  Be(m)  • 

If  -  I  fl  <  (2w  )_2C  ( 1tCd)”**<n^  If  “I  +  E  IL*(*+2wj)l  ,.lf*l. 

n  •  1V  '  1  jSO  "  »,supp  f  1 

<  [c  (1+Cd)’M(n)  +  Z  <1+C|j|)'M<n)]  If “I 

je  \  oujuj’ 

<  C1(1+Cd)'M(n,lf*l1  .  Ml 


We  now  discuss  briefly  the  particularly  symmmetric  and  special  case 


n  “  (r,r,r)  . 


Figure  5.2  shows  the  level  lines  for  P_  :*  p.  _  .  for  r  =  3  .  Note  that  its  minimum 

r  l  *  #  *  #  *  J 

seems  to  occurs  at  (2w/3,2*/3)  ,  and  this  can  be  verified  analytically  for  r  <  4  .  We 
conjecture  that  this  is  no  accident,  but  is  the  case  for  all  r  . 

Figure  5.3  shows  L_  L,_  _  for  r  -  3  (?)  .  The  fast  decay  is  ouite  striking, 
making  plain  that  cardinal  interpolation  with  this  L  would  be  strongly  essentially  local. 

If  we  assume  that  f*  e  Lj  with  supp  f*  c  S)^  ,  then  the  convergence  of  Ipf  can  be 
stated  in  a  particularly  nice  way.  He  define  a  bivariate  "Whittaker"  operator 

Wii2— >  L2sf  |— >  r  2  f(j)x*<— j)  , 


with 


a  cos  -§^(u+v)  COS  -|^(2u“v>  COB  t^(v-2u) 
X*(u«v)  [  - 3- -  ♦  - 2 -  ♦  - 2 -  ] 


2« 


(u-2v)(v-2u)  (u+v) (u-2v) 


(u+v) (v-2u) 


the  Fourier  transform  of  x  s*  X|t  (  •  Note  that  the  translates  of  x*  are  orthogonal 

in  Kj  .  As  in  the  univariate  case  [S21  ,  the  "Whittaker”  series  provides  the  limiting 

operator  for  cardinal  interpolation  I(rtrtrj  as  r**»  .  More  precisely,  we  have 

1W  -  I  :  l, — >  LI  - >  0  .  (5.7) 

r  2  2  r+~ 

It  t  e  t~  and  supp  f"  CO.,  ,  ..  ,  we  have  w((f(j))  ,)  »  f  and  hence  (7)  is  an 
Lj -version  of  Theorem  5.2. 

To  proem  (7),  we  first  show  that  the  cardinal  interpolation  maps  I_  :»  I.  _  .  are 

bounded  as  maps  from  &2  to  Lj  uniformly  in  r  :  For  f  e  l2  ,  we  have 

IIrfl2  -  IIf(j)e"ijVl2  <  lfl2(l  lL*(.+2»j)l^n  )  <  C  lfl2  . 

By  the  uniform  boundedness  principle,  it  is  therefore  sufficient  to  check  the  convergence 
(7)  for  the  unit  vector  Cj  €  i2  .  But  this  is  an  obvious  consequence  of  (3).  ||| 


In  the  univariate  setting  [R]  proved  the  convergence  in  1^  for  all  p  e  ( 1 , <•)  ,  but  we 
leave  the  corresponding  bivariate  problem  to  a  later  paper. 

[MRR]  extend  Schoenberg's  univariate  result  to  include  the  possibility  that  f  has 
support  at  ±».  Wiis  requires  the  realization  that  Lr<±x)  — >  1/2.  Theorem  5.1  says 


-  5.7 


nothing  about  the  limit  of  L“(x)  in  case  x  e  3ft  .  For  the  special  choice  n  -  (r,r,r) 

r  m 

such  a  statement  is  relatively  easy  to  make. 


Corollary 


For  the  special  choice  n 


lim  L*(x) 
n 


(r,r,r )  , 

1/3  ,  x  e  Az 
1/2  ,  x  e  38  \  Az 


with  z  (2*/3,2*/3) 


(5. 


Clearly,  our  result  concerning  the  convergence  of  In  still  holds  if  f  is  a 

measure,  absolutely  continuous  in  a  neighborhood  of  38  and  supported  in  8  .  Our 

m  m 

result  is  best  possible  in  the  following  sense.  If  supp  f*n3n  f  0  ,  then,  in  general, 

Inf  does  not  converge  to  f  .  For  example,  if  f(x)  *>  cos(z*x)  ,  then 

I  f(x)  - >  I  e****'  x  =  [  cos  ■|~(u+v)  +  cos  t~(2u-v)  +  cos  -§^(2v-u)  ]  . 

r  r*«  AeA  1  3  3  3  ‘ 

This  follows  from  the  Corollary.  However,  it  requires  slightly  more  precise  information 
about  the  convergence  asserted  in  (8).  The  heuristic  argixnent  is  that 

"  V  Ln  *  1  (2,°2  E  «.*2.jt;(-2W3» 

and  therefore 

W  ,/j)  • 


Here,  6  denotes  the  Dirac  measure  at  £ 


6.  Detailed  eetiaatee.  In  this  section,  we  prove  Propositions  5.1  -  5.4  and  various 

lemmas  needed  in  the  proofs,  this  amounts  to  a  detailed  study  of  the  functions  an  j  and 

the  set  0  and  how  they  depend  on  n  @  N  .  In  particular,  we  need  to  study  the  boundary 
n 

of  .  This  boundary  is  made  up  of  pieces  of  curves  given  implicitly  by  the  equation 


a  j(x)  «  1 
n,  1 


for  some  )  8  J 


We  use  the  symmetries  of  the  given  situation.  Recall  the  notation 


o(n)  i-  n0(2)'  n<,(3)>  * 


We  conclude  from  (3.5)  or  directly  from  (2)  that 


an,3(x)  *  a«(n>,A;j(V>  • 


This  implies  that 


and  therefore 


where,  to  recall  from  Section  4, 


»*fo  np)  «  (l  n  r  , 

o'  n  '  o(n)  o 


Next,  we  conaider  the  boundary  of  Q  (cf.  Figure  5.1).  Each  -)  e  J  lies  in  two 

n 

cone*  R  •  With  C_j  their  union,  we  define  the  curve 
<J  j 

r  ,  :■  {x  e  C  a  ,(x)  »  1}  .  (6 

n,  3  -j  n,  j 

The  boundary  of  is  made  up  of  segments  of  these  curves.  It  follows  from  (1)  that 


A  r  . 
a  n,  j 


o(n),A*j 


6.1.  For  j  e  J  ,  denote  by  )0,  the  vectors  spanning  the  union  Cj  of  the 


two  cones  Rfl  containing  j  .  For  n  e  (0,*»)  ,  the  curve  passes  through  the 

points  -jg,  -j/2,  -jj  and  is  syiwetrlc  with  respect  to  the  point  -j/2  .  Moreover,  it  is 
monotone  (in  an  appropriate  coordinate  system)  and  is  Lipschits  continuous,  uniformly  in 


6.1 


Proof.  The  symmetry  with  respect  to  the  point  -j/2  follows  from  the  relations 


an,j(x-3>  “  1/an,-j(x>  =  1/an,j'-x) 

which  are  immediate  consequences  of  the  definition  (5.2)  of  an  .  For  the  rest, 
view  of  (5),  sufficient  to  consider  j  =  (-1,0)  .  In  this  case,  the  curve  T 

n,  g 

by  the  equation 


Using  the  fact  that  u  >  0 


|_“_)r|_i!±v_|t  =  , 

1 u- 1  1 u+v- 1 1 


u+v  >  0  ,  and  solving  for  v 


we  obtain 


it  is,  in 
is  given 


(6.6) 


v  =  -u  +  1/[l  ♦  (-^)°]  ,  0  <  u  <  1  ,  (6.7) 

where  a  :=  r/t  .  This  shows  that,  for  any  q  e  (0,®)  ,  the  points  -jg  =  (1,-1),  -j/2  = 

(1/2,0),  and  -j1  ®  (0,1)  lie  on  the  curve.  Moreover  we  have 

dv/du  -  -1  -  [l  ♦  M-u)'2  (6.8) 

which  shows  that 


dv/du  <  -1  ,  (6.9) 

with  equality  only  if  u  -  0  .  this  proves  the  remaining  assertions  of  the  Lemma.  | | | 


6.2.  Let  jQ,  j.,  be  the  two  vectors  which  span  the  cone  R  .  For  n  e  (0,®)  , 


the  curves  T 


and  T  intersect  at  a  unique  point  z  e  R  .  The  boundary  of 

*J0  -  n,-j1  - * - -  n,o  o  - “ - 


fln  consists  of  the  segments  of  the  curves  Tn  ^  connecting  the  points  in  J/2  and 
n,o  3 


Proof.  In  view  of  (3)  and  (5),  it  is  sufficient  to  consider  the  case  a  ®  (1),  R  = 

o 

R  and  jg  =  (1,0),  =  (0,1)  .  By  (9),  Cn  (  1  qj  has  slope  <  -1  ,  with  equality  only 

at  the  point  (0,1)  .  Similarly,  a  direct  computation  shows  that  T  .  .  has  slope 

n, (0,-1) 

between  -1  and  0  .  Since  T  ,  ,  and  r  ,  „  ..  pass  through  the  points  (0,1), 

n,  (-1,0)  n, ( 0,-1 ) 

(1/2,0)  and  (0,1/2),  (1,0)  respectively,  they  intersect  at  a  unique  point  z  e 

n,  (  1 ) 

2 

(0,1/2)  .  To  show  that  the  boundary  of  0^  intersected  with  R  consists  of  the  segments 

connecting  (0,1/2)  with  zn^ ( 1 )  •  and  zn,(1)  with  (1/2,0)  ,  we  prove  that,  for  x  = 
(u,v)  e  R  , 

an  (-1  0)^x^'  an  (0  -1)*x^  <  1  implies  that  an  ^ ( x )  <  1  for  all  j  e  J  . 


2 


Indeed, 


■  '^T'r'^rr't  <  1 

implies  u  <  1/2  ,  and,  from  an#(o,-i)(x)  <  1  .  It  follows  that  v  <  1/2  .  This  implies 

V<-i,i)(x)  "  'l^T1  ,^+Tl  <  1  ' 

and  the  other  cases  can  be  checked  just  as  easily.  III 


6.3.  For  all  m  e  N  ,  T  ,  -  Hm  {T  n  e  (0,-)3ntl}  .  Hence  Lemma  1  is 

-  m,  j  n-nn  3 


valid  for  all  n  e  N  . 


Proof.  Without  loss  of  generality,  we  consider  only  the  case  j  »  (-1,0)  .  We  claim 


BL(  (0,1),  (0,1/2),  (1/2, 0) ,  { 1 , -1/2),  (1,-u)  if  m^“  0  ai 
fl  BL((0,1),  (1/2, 1/2),  (1/2,0),  (1/2, -1/2),  (1,-1))  if  m^  - 


and/or  m 


and/or  m3«  0 


where  BL(Xj, . . .  ,xm)  denotes  the  broken  line  with  vertices  x^,  xR)  ,  Consider,  e.g. 

the  first  case.  By  (9)  and  the  symmetry  of  the  curves  T  .  and  T  .  ,  we  have  for 

n, j  m, j 


n  e  ( 0 ,“)  that 


dist(T  ,  r  )  <  u  ,  with  u  s.t.  a  (u  ,1/2)  -  1 

m,3  n,3  n  n  nr  j  n 


From  (6),  we  obtain 


ln(1-un>  -  In  un 

T  S“  t/r  "  ln(u  +1/2)  -  ln( 1/2-u  )  * 

n  n 


If  n  — >  m.  with  m^  »  0  and/or  m3  »  •  ,  we  must  have  y  — >  •  and,  by 


the  above 


equation,  this  implies  that  u  — >  0  .  The  second  case  can  be  handled  similarly.  I  I  I 

n 


Lemma  3  is  a  particular  case  of  the  next  lemma  which  states  that  T  ,  depends 

n«  J 


continuously  on  n  . 


6.4.  For  n,  m  e  N  ,  dist<r  ,  ,T  ,)  - >  0  as  n  — >  m  . 

-  n,j  m,  j  — 


Proof.  In  view  of  Lemma  3,  we  may  assume  that  m  e  (0,**)  .  Moreover,  it  is  sufficient 


to  consider  j  ”  (-1,0)  .  In  this  case,  it  follows  from  (7)  that  T  .  — >  T  .  pointwise, 

n,J  m,j 


-  6.3 


both  curves  being  viewed  ss  functions  of  u  •  By  the  uniform  Lipschitz  continuity  of  the 
curves,  this  implies  the  assertion  of  the  lemma.  I | | 

Lemmas  1-4  give  a  qualitative  description  of  the  boundary  of  .  We  summarize  the 
main  features  in  the  following 

Corollary.  (1)  3fln  consists  of  segments  of  the  curves  Tn  connecting  the  points 

in  J/2  with  the  intersection  points  z  #  o  e  t.  • 

-  —  -■  ■  -■  - —  n,o  3 

(11)  3(1  Is  piecewise  monotone  and  is  Lipschltz  continuous,  uniformly  in  n  . 

(ill)  depends  continuously  on  n  in  the  Hausdorff  topology. 

Note  that  this  provides  the  proof  of  Proposition  5.1. 

To  give  a  few  examples,  we  list  below  all  cases  for  which  $1^  has  a  piecewise  linear 


boundary 


Moreover,  we  have 


(1,1,0) 

*(1,0,1) 


(-,1,0) 

’(-.O.l) 


(1,-,0) 

Jd,0,-) 


(6.10) 


(0,1,1) 


<0,-,1) 


(0,1,-) 


We  take  the  occasion  to  prove  the  following  observation  which  stresses  the  underlying 
hexagonal  structure. 


Proposition  6.  $).  :»  Ofl  -  int  conv  J/2  ,  £),  LJfl  “  A  (0 , 1/2 ) 2  . 

On  i  n 


(1/2, 1/2) 


Figure  6.4 


Proof.  We  claim  that,  for  any  n  , 


{(u,v)  s  u,  v  >  0,  u+v  <  1/2}  C  R  n  JJ  c  [0,1/2r  . 

"  n  “ 


(6.11) 


This  follows  from  Lemma  1,  in  particular  from  the  fact  that  the  curves  T  ....  and 

n,  (-1 , 0) 

r  pass  through  the  points  (1/2,0)  and  (0,1/2)  and,  as  functions  of  u  ,  have 

n,  ( u,-l j 

slopes  <  -1  and  »  -1  ,  respectively.  To  complete  the  proof,  note  that  n  »  (0,1,1) 


6.6 


gives  equality  in  the  first  inclusion  of  (11)  while  n  =  (1,1,0)  gives  equality  in  the 
second.  | | | 

We  are  now  also  ready  for  the 

Proof  of  Proposition  5.3.  Because  of  the  continuity  of  0^  as  a  function  of  n  ,  it 

is  sufficient  to  consider  n  e  (0,«)3  .  In  this  situation.  Figure  5.1  gives  a  qualitatively 

correct  description  of  Si  .  Because  of  the  geometry  of  SI  and  the  symmetry  relations 

n  n 

(2),  it  is  sufficient  to  establish  the  following  claims: 

<i>  j  ♦  r  .  -  r  .  for  all  j  e  J  . 
n,  3  n»"j 

(ii)  The  curve  (1,0)  +  T  ,  passes  through  the  point  z„  ...  . 

n, ( 1,-1 )  n, '  1  < 

The  first  assertion  follows  from  the  relation 

an,j<x'j)  *  1/an,-j(x> 

alluded  to  earlier  and  directly  derivable  from  the  definition  (5.2)  of  an,j  •  A8  to  (ii), 
note  that 

an,(-1,0)(*n,(1))  “  1  “  an,(0,-1)(zn,(1>> 

implies  that 

1  "  an,(0,-1)/an,(-1,0)<zn,(1)) 

“  an,(1,-1)(zn,(1)  "  (1'0,)  ' 

zn, ( 1 )  e  <1'0)  +  rn, ( 1,-1 )  '  1,1 

The  next  three  lemmas  state  various  estimates  for  the  functions  a„  .  needed  for  the 

J 

proof  of  Proposition  5.2. 


Idmtmm  6.5.  For  n  6  [0,*»)  with  at  most  one  component  less  than  1  ,  we  have 

a  <  [i  +  c  dist( x,  r  .ns  )1  1  ,  x  e  fl  ,  )  e  j  , 

1  n,j  n  1  n 

with  c  a  positive  constant  which  does  not  depend  on  x  ,  n  ,  o£  j  . 


Proof,  we  may  assume  tnat  x  ««  (u,v)  e  R  ,  in  particular  that  u,  v  e  [0,1/2]  .  We 

consider  each  j  e  J  separately  and  suppress  all  references  to  n  . 

r  t 

(i)  j  «  (1,0)  .  We  have  a^fx)  ”  I  ^ I  <  1/2  .  This  proves  (12)  since 


dist(x,  T  .n  a  )  <  1  .  For  the  estimate,  we  have  used  the  fact  that  min  {r,t}  >  1  • 
n,  J  n 

The  case  j  «  (0,1)  is  similar. 

(ii)  j  *»  (-1,0)  .  The  following  figure  may  be  of  help  in  following  the  argument. 


Figure  6.5 

Let  i  ;«  (0,-1)  and  let  z  »:  (u0,v0)  :«  be  the  intersection  of  the  two  curves  T 

and  r  .  We  consider  two  cases. 

For  v  <  vQ  ,  there  exists  e  >  dist^lx,  1*^  n  fl)  such  that  (u+e,v)  e  I\  ,  i.e., 

aj(u+e,v)  ■  1  .  It  follows  that 

a^(u,v)  -  a^ (u,v)/a j (u+e,v) 

/  u  ■ir /1-u-e-ir r  u+v  \t; 1-u-v-e\t 
'•u+e'  '  1-u  '  '  U+ V*  6  '  '  1  -u-v  ' 

<  f_2lv_)«wx{r,t}  <  +  -1  _ 

'll  +  V-f  c  * 


where  we  have  used  the  fact  that  s,  t  >  1  and  the  last  inequality  is  easily  checked. 

Since 

|u-u0l  <  I u |  +  |uQ|  <  2 | u |  , 

this  proves  (IS)  for  this  case. 

(b)  x  e  r  .nil  *  l-e.#  0  <  u  <  u,  =  |u.|  .  We  assume  first  that 

vQ  -  1/2  >  lu0l/3  .  (6.18) 

Since  z„  e  Tj  ,  we  have 


a^(x) 


* j  <  *o)^a  j(x)  “ 


> 


.  ( U+ 1 ) /u  . r ,  (1-v)/v 

(V1)/uo  (1“V0)/V0 


I  1 


♦ 


(V0~V>/<VV0> 

|1-vo)/vo 


I 


s 


> 


+  (v0-v)  . 


In  view  of 


|u  -  uQ|  <  2|uQ|  <  6(vq-1/2)  <  6(vq-v)  , 
this  proves  (15)  under  the  assumption  (18). 

Next,  suppose  that 

vQ  -  1/2  <  |uq|/3  .  (6.19) 

We  claim  that 

|(u0+1)/u0|r  >  1  +  c | Up  I  ,  (6.20) 

and  this  finishes  the  proof  of  (15)  for  this  case,  in  view  of  (16)  and  the  inequalities 

1  r  u  +1  r 

1/a^(x)  >.  I^^-I  >  |— — |  >  1  +  c|uQ|  >  1  +  (c/2)  |u  -  uQ I  . 

To  prove  (20),  we  use  the  fact  that  z0  ft  Q)  ”  3  »  therefore 

1  ( u0+ 1 )  /u0  T 1  ( Vv0 )  /  ( 1  -  V  vQ ) 1 1  “  1  • 

Solving  this  for  r  ,  we  obtain 


r 


In  l<«0»v0>/M-u0-v0)l 

t  In  |(u0+1)/u0l  * 


1/2  -  (2/3) I uQ I | 
1/2  +  (2/3) |uQ| I  / 


In  |(u0+1)/u0l  . 


Here  we  used  the  assumption  (19)  and  that  t  >  1  .  Therefore,  we  have 


r  In  |(u0+1)/u0l 


1  +  (4/3) |u0 | 
1  -  (4/3 ) | u Q | 


>  In  |1  +  (4/3) | Uq | | 


6.11 


which  establishes  (20).  This  completes  the  case  j  *  (1,-1) 
The  case  j  »  (-1,1)  is  treated  similarly.  ||| 


Proa  the  statement  of  Proposition  S.2,  recall  the  definition 

J'  i-  A*(1,1)  -  {*(1,1),  *(2,-1),  *(-1,2))  . 


6.6.  For  n  6  (0,»)  on,  x  e  8  ,  and  j  6  J' 

n  1 


i-i 


•n,j<xl  *  f1  *  c  dlsttx,  -J/2 ) ]" 
with  C  a  positive  constant  which  does  not  depend  on  n  ,  x  ,  oz_  j 


'l-u'  'l-v'  2-u-v  1/2+e  'l+t' 


Proof.  Assume  without  loss  that  x  -  (u,v)  6  R  . 

We  consider  first  J  «  (-1,-1)  .  Assume,  e.g.,  that  1/2  -  u  -  dist^lx,  (1,1)/2) 

•i  C  .  Then  we  have 

*(-1,-1),x) 

In  the  reaulnlng  cases  j  e  J'  \  (-1,-1)  ,  we  have  dlst^l-J^.R)  <  2  and  therefore  it 
is  sufficient  to  bound  a^tx)  by  an  absolute  constant  less  than  1  .  This  is 
straightforward,  using  the  fact  that  »(.1i0}(*>  ,  *(0,-i)(x>  list  only  the 

estimates i 

i  m  <V>  » 


r  s  t 

.  u  .  ,  v  ,  i  u+v  ,  ,  ,-r ,-s  -t 

*u+1*  v+ 1  u+v+2  <332 


1  ■  (2,-D  > 


*(2,-1)(x)  *  a(2,-1)(x>/*(0,-1)(x)  "  '2+u1  ^ 2-u-v ^  <  5  2 


3  '  <-2>D  » 


•(-1,2)(X)  *  “(-1,2)(x,/*(-1,0)(x) 


2tv  2-u-v 


The  remaining  two  cases  ,  j  »  (-2,1)  ,  (1,-2)  are  similar.  Ill 


6.7.  For  n  «  N  ,  x  e  fl"  ,  and  1  «  S2  \  (OUJUJ’)  , 
■  “  n  ■■  ■ 


•n> j(x)  <  [l  ♦  cl) | l”1  , 


with  C  a  positive  constant  which  does  not  depend  on  n  ,  x  ,  or  J  . 


Proof.  Let  j  »  (k,t)  and  x  »  (u,v)  and,  without  lose  of  generality,  x  e  R  , 
hence  u,  v  e  [0,1/2]  .  We  consider  several  cases. 

(i)  k,  t  f*  0,  -1  ,  and  k+t  ^  0,  -1,  -2  i  Then 


1^1  <  (i  ♦  Mel)-1  ,  l^fl  <  n  ♦  I*!)"1  .  Iu+v+k+t 


U+V  I  <  <1  ♦  lk+tl/3)-1 


Using  the  fact  that  ( 1+p)-1 ( 1+q)”1  <  (1+p+q)”1  for  p,q  >  0  and  that  |p|  ♦  |p+q| 

>  max  {|pi,lq|}  ,  we  see  that  the  product  of  any  two  of  the  above  lefthand  terms  is 
bounded  by  (1  +  ljl^/3)"1  .  Since  at  most  one  of  r  ,  s,  t  is  less  than  1  ,  this  yields 


an,j(x)  *  <1+Ht|)-r<1+lt|)‘S(1  +  |k+t|/3)-t  <  (1  +  Ijly3)_1 


(11)  j  “  (0,-2)  >  Since  x  (  ,  we  have 

,<*>  -  l^|8|,r^|t  <  1 


“n, (0,-1 ) ' 


and  therefore 


■  1-V.B.  l-u-v,* 


an,j<*>  <  *n,j<X,/an,<0.-1)<>0  -  1^1  1^1  <  2"82_t  <  1/2  -  (1  4  Ij^)’1 


(ill)  k  -  0  ,  l  *  -2, -1,0  :  Here 


v  u+v  t 


«„,i(x)  ■  '^i1  <  d+uD^d+ui/s)*1  <  n  ♦  uiys)"1 


(iv)  k  -  -1  ,  l  *  Op  1f  2  .  Here 
r 


■  u+v  fc 


-1 


1^1  ^  *  ’•<iHt|)-,(i+it-i|/3)-t  <  (1  +  Ijiys) 

(v)  k  -  0,  -1  and  1  “  0,  -1  s  Treated  in  analogy  to  cases  (ii)  -  (iv)  . 


(vl)  l  ■  -k  ,  k  +  -1,  0,  1  x  Here 


•»  •!<*>  -  1^:1  hrd  <  (i+iki)-r(i+iki)'8  <  d  +  ui j"1 


n, J'“*  'u+k'  ’v-k 

(vii)  t  -  -k-1  ,  k)i-2,  -1,  0,  1  s  Here 


,1-u,r.  v  8 


8n,j<*>  <  an,J<x>/an,(-1,0)<*>  <  ^  '^T1  «  |k|  (1+|k-,|)  <  d+IM/2) 

<  (i  +  ijiy3)_1  . 


-i 


(viii)  l  «  -k-2  ,  k  f  -2, -1,0  :  Here 


.  r  s  .  .  t 

u-1 ,  , v  .  iU+V-1 , 


-1 


an,j<*>  4  'n,j<xl/in,M,0),xl  “  W  fe1  4  (1+|k|/3) 


<  (i  ♦  uiy9)‘ 


This  covers  all  j  e  S2  \  (0UJUJ' )  .  ||| 


Proof  of  Proposition  5.2.  Lemmas  5-7  prove  (5.4)  for  n  e  N  n[0,»)  ,  with  the 


constants  independent  of  n  .  Since  I)  depends  continuously  on  n  e  N  (by  Proposition 


msmavsgK 


5.1),  this  proves  (5.4)  for  all  n  e  N  .  Ill 


Proof  of  Proposition  S.4.  To  prove  (5.6),  we  consider  three  cases. 

(i)  j  e  J  :  From  the  geometry  of  the  set  (1^  =  (ln,  (cf.  Figure  5.1  which  gives  a 

qualitatively  accurate  description  of  the  general  situation)  and,  in  particular,  from  the 

estimates  of  the  slopes  of  the  curves  T  .  ,  we  can  see  that  in  this  case 

n,j 

dist(x,3Sl  )  =  dlst(x',r  .  n  (1  )  . 

n  n ,  ]  n 

Therefore,  (5.6)  is  a  consequence  of  Lemma  6.5  in  this  case. 

(ii)  j  e  J*  :  Assume,  e.g.,  that  x'  e  Rnil^  .  Then,  for  j  e  J'  \  (-1,-1)  ,  we  have 
dist— (x' ,-j/2)  >  1/2  and  (5.6)  follows  from  Lemma  6.6  since  dist^tx, 33  )  <  C  . 

It  remains  to  consider  the  case  j  =  (-1,-1)  .  From  the  bounds  on  the  slopes  of  the 

curves  V  ,  „  T  »  we  see  that  Q  lies  in  the  half  space 

n, ( 1 ,0)  n, (0, 1 )  n 

{y  s  y* ( 1 , 1 )  >  -(u0.v/(1,1)} 

where  (uq,Vq)  s»  is  the  point  of  Intersection  of  the  curves  r n  and 

r  ...  Since  (1  flB  C  [0,1/2]2  ,  it  follows  that,  for  x‘  e  0  OR  , 

n,  ( 0 ,  •  1 )  n  n 

diet,  (x, 3(1  )  -  dist,(x',  -*.,.+(1,1)1  .  (6.21) 

i  n  1  (l) 

Here,  dist  ,  denotes  the  t ,-distance.  Moreover,  we  have,  in  view  of  l“Z(i)  +  (1#D)  ~ 

(1/2, 1/2)  -  (1/2, 1/2)  -  z(1)  ,  that 

dist  (x',  -*.,,+(1,1))  =  dist ,( x' ,z, , , )  +  2dist .(z..,, (1/2, 1/2)) 

1  (1)  1  (1)  1  (1)  (6.22) 

<  2dist,(x\-j/2)  . 


This,  together  with  (21)  and  Lemma  6,  proves  (5.6)  for  this  case. 


(iii)  j  e  X2  \  (OUJUJ' )  s  In  this  case,  (5.6)  follows  from  Lemma  7  since,  for  any 


j  *  0  , 


disttx,  3(1  )  <  C |  j  |  . 

n 


This  completes  the  proof  of  Proposition  5.4.  | | | 
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